Measurement of anisotropic radial flow in relativistic heavy ion collisions 
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We suggest the azimuthal distribution of mean transverse (radial) rapidity of the final state par- 
ticles as a more direct measure of the transverse motion of the source than the standard azimuthal 
multiplicity distribution. Using a sample generated by the AMPT model with string melting, we 
demonstrate that the azimuthal amplitude of the suggested distribution characterizes the anisotropic 
radial flow, and coincides with the parameter of anisotropic radial rapidity extracted from a gener- 
alized blast-wave parametrization. 
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I. INTRODUCTION 

Relativistic heavy ion collisions provide a vi^ay to study 
the properties of strongly interacting matter. The ob- 
servation of large elliptic flow at RHIC is considered 
one of the most important signatures for the forma- 
tion of the strongly interacting Quark Gluon Plasma 
(sQGP) P, m. The Row harmonics are Fourier coeffi- 
cients of the azimuthal multiplicity distribution of final 
state hadrons Q . They are considered sensitive probes of 
the evolution of the system formed in relativistic heavy 
ion collisions [3]. 

One common feature of flow harmonics is their mass 
ordering in the low transverse momentum region 
This phenomena can be well understood by hydrody- 
namics with a set of kinetic freeze-out constraints, i.e., 
the temperature, the radial flow, and the source defor- 
mation Q- The radial flow is usually described by 2 
parameters. The first is the isotropic radial velocity (or 
rapidity, vt ~ tanhyx)- It presents the surface profile 
of isotropic transverse expansion of the source at kinetic 
freeze-out. 

The other parameter is the anisotropic radial velocity 
(i.e., the azimuthal dependent radial velocity). It mea- 
sures the difference of the radial flow strength in and 
out of the reaction plane. It is introduced to account 
for the anisotropic radial flow field which arises in non- 
central collisions. The observed elliptic flow can be gen- 
erated by anisotropic radial flow Q. Moreover, the 
shear tension of viscous in hydrodynamics is supposed to 
be proportional to the gradient of radial velocity along 
the azimuthal direction Q, which is directly related to 
anisotropic radial velocity. The proportionality constant 
is the shear viscosity. 

In hydrodynamic model [8l-ITo|. these parameters are 
not independent. They are related by the initial condi- 
tions and the equation of state and their determination 
is crucial for theoretical calculations. 

The azimuthal distribution of the mean transverse ra- 
pidity of final state hadrons ((yT(0))) directly measures 
the transverse motion of the source at kinetic freeze- 
out [ll| . It should be helpful in determining the parame- 



ters of the anisotropic radial rapidity. In contrast to the 
azimuthal multiplicity distribution, where only the num- 
ber of particles is concerned, here the average is over all 
particles in a given azimuthal direction. The influence of 
the number of particles is excluded. 

As we know, {yT^cj))) should contain three parts: av- 
erage isotropic radial rapidity, average anisotropic ra- 
dial rapidity, and average thermal motion rapidity Q. 
Since both thermal and radial motions contribute to the 
isotropic rapidity of the distribution, the isotropic ra- 
dial rapidity itself can not be directly obtained from 
the distribution. So conventionally, the radial flow 
parameters are extracted from the pj spectra of the 
hadrons [H, [ll], or dileptons [l^, by generalized 
blast- wave parametrization 0, Q. 

Fortunately, the thermal motion is isotropic. As such, 
it docs not contribute to the anisotropic radial flow. The 
azimuthal amplitude of the mean transverse rapidity dis- 
tribution should correspond directly to the anisotropic 
radial rapidity. It is interesting to see the features of the 
azimuthal distribution of mean transverse rapidity, and 
how its azimuthal amplitude relates to the parameters 
of anisotropic radial rapidity extracted by a generalized 
blast-wave parametrization. 

In the paper, we define the (yT(<?^)) in section II. Us- 
ing a sample generated by the AMPT model with string 
melting [T^ , the suggested distribution and its parti- 
cle mass and ccntrality dependence are presented. These 
show that the isotropic and anisotropic parts of the sug- 
gested distribution behave as the expected radial flow 
(with a random thermal component), and anisotropic ra- 
dial flow, respectively. In section III, the spectra of 
6 particle species and their corresponding elliptic flows, 
V2{pt), are presented. Fitting these spectra and ellip- 
tic flows by a generalized blast-wave parametrization, 
the temperature, and the radial flow parameters are ob- 
tained. It is found that the parameter of anisotropic 
radial rapidity is well described by the azimuthal ampli- 
tude of the suggested distribution. Finally, the summary 
and conclusions are given in section IV. 
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II. AZIMUTHAL DISTRIBUTION OF MEAN 
TRANSVERSE RAPIDITY 

Usually, the transverse rapidity of a final state hadron 
is considered a good approximation of its transverse ra- 
pidity at kinetic freeze-out It is defined as. 



j/T = ln( 



rriT +PT -. 
Too 



(1) 



where toq is the particle mass in the rest frame, 
is transverse momentum, and tot — \/ TOq + y\ is the 
transverse mass. The mean transverse rapidity in a given 
azimuthal angle bin is defined as the summation of all 
particles' rapidities divided by the total number of par- 
ticles, i.e., 
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yT,i(</'m-V'r), (2) 



where t/|n ^ is the transverse rapidity of the ith particle 
and iV^ is the total number of particles in TOth azimuthal 
angle bin in eth event. The direction of reaction plane 
is given by V'r, which is zero in model calculation, and 
can be determined in experiment by 3 similar ways as 
that for the azimuthal multiplicity distribution [? ] . By 
definition, the influence from the number of particles is 
therefore removed. Eq. (2) measures the mean transverse 
motion in azimuthal direction 

In order to see the generic features of the defined dis- 
tribution, we present the distribution from a sample of 
Au -I- Au collisions at 200 GeV generated by the AMPT 
model with string melting in Fig. 1(a). For the fol- 
lowing comparison, we use a longitudinal rapidity win- 
dow of Ij/lI < 0-1: the same as published spectra data 
from the STAR experiment [ll,!^]. 

From the figure, we can see it is a periodic function of 
azimuthal angle and can be well fitted by. 



(yT(0)) = 2/TO + yT2 COS(20). 



(3) 



Eq. (3) consists of two parts: an isotropic mean rapidity, 
2/To = 1.3371 ± 0.0001, and a mean azimuthal dependent 
rapidity amplitude, y^i = 0.0334 ± 0.0002. This directly 
indicates the anisotropic distribution of transverse mo- 
tion, in addition to the known anisotropy distribution 
of particle number. Further, the anisotropic amplitude, 
?/T2j should correspond to the parameter of anisotropic 
radial rapidity. 

The isotropic part is a combination of radial rapidity 
and thermal motion rapidity. As we know, the thermal 
motion is mainly determined by the temperature and 
particle mass. For a system at fixed temperature, the 
lighter particles should have larger thermal velocity. To 
see this feature in isotropic rapidity we plot the distri- 
butions of three different particle species, tt, fc, and p 
in Fig. 1(b). It shows that the lightest particles (pions) 
have the largest isotropic rapidities, while the heaviest 



particles (protons) have the smallest isotropic rapidities, 
and intermediate mass particles (kaons) have rapidities 
between them. These indicate that their isotropic ra- 
pidities are ordered as expected from random thermal 
motion. 

To see the infiuence of centrality on the anisotropic 
part, the azimuthal distributions of mean transverse ra- 
pidity at three typical centralities, 0—5%, 30%— 40%, and 
60% — 70%, are presented in Fig 1(c). From the figure, we 
can see that the distributions are azimuthal angle depen- 
dent in non-central collisions, i.e., the mid-central and 
peripheral collisions with centralities of 30% — 40% and 
60% — 70%. It becomes almost fiat and azimuthal angle 
independent in central colhsions (0 — 5%). So the az- 
imuthal dependent part of the distribution appears only 
in non-central collisions. 

This is consistent with the fact that only two parame- 
ters, the temperature and radial rapidity, are required to 
describe the observed px spectra in central collisions, as 
done in early blast- wave paramctrization Q. However, 
the parameter of anisotropic radial flow is necessary for 
non-central collisions 0, Q . 

The centrality dependence of the azimuthal amplitude, 
2/T2, is shown in Fig. 1(d) by red solid circles. It has 
a maximum in mid-central collisions, decreases toward 
peripheral and central collisions, and is close to zero in 
central collisions. 

The disappearance of ?/t2 in central collisions also in- 
dicates that the thermal motion, which exists in central 
collisions as well, is isotropic and does not contribute to 
the anisotropic radial rapidity. Therefore, the azimuthal 
amplitude of the suggested distribution describes the pa- 
rameter of anisotropic radial rapidity. In order to show 
this quantitatively, we will compare it with the parameter 
that is extracted from the same sample by a generalized 
blast-wave paramctrization in the following section. 

III. THE PARAMETERS OF RADIAL FLOW 

The blast-wave model is currently the only model that 
simply includes the radial flow parameters. It is mo- 
tivated from hydr ody namics with the kinetic freeze-out 
parameters d, [l3, E, [2l| - [23j . It is assumed that the lon- 



gitudinal expansion is boost invariant j24l |. The single- 
particle spectrum is given by the Cooper-Frye formalism 
(as in hydrodynamics) j25| . 
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p^daf,{x)f{x,p), 



(4) 



where f{x,p) is the momentum distribution at space- 
time point X. Eq. (4) is an integral over a freeze-out 
hyper-surface, and sums over the contributions from all 
space-time points. 

Originally, local thermal equilibrium is assumed to be 
reached at kinetic freeze-out and a Boltzman distribution 
of the momentum is applied [l2j. It has been shown re- 
cently that a Tsallis distribution provides an even better 
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FIG. 1: (Color online) (a) {yT{4>)) of minimum bias sample, {h){yT{4>)} for 3 different mass particles, (c) (jjT^ij))) at three 
different centralities, — 5% (red points), 30% — 40% (black triangles), and 60% — 70% (green down triangles), (d) The 
centrality dependence of yT2 (red points), and p2 (black triangles). 



description for all spectra from elementary to nuclear 
collisions [l^ll^l- So, we use the Tsallis distribution for 
f{x,p), i.e.. 



g-1 
T{x) 



p ■ u{x) — fi{x) 



(5) 



where q is the parameter characterizing the degree of non- 
equilibrium, and T is the kinetic freeze-out temperature. 
Thus the transverse momentum spectrum can be given 
by m, 



source element defined with respect to the reaction plane. 
They are related by tan(0(,) = (i?jf/i?y)^ tan((/)s). 

There are 5 undetermined parameters: the tem- 
perature (T), isotropic radial flow rapidity (po) and 
anisotropic radial flow rapidity {p2), 9 of the Tsallis 
distribution, and Rx/Ry- Since all the particles are 
assumed to move with a common radial flow velocity, 
the mean kinetic freeze-out parameters are usually ob- 
tained by the simultaneous fitting of spectra from sev- 
eral hadrons [13, [U and elliptic flow [6| . Elliptic flow, 
'V2{pt), is the second coefficient of the Fourier expansion 
of azimuthal multiplicity distribution (32l . [ssj . and de- 
fined as. 



dN 



PTdpi 



/ d(l)s / dye^^b ^ coshj/ / mrrdr 
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-PT sinh p cos((/)f, — 



(6) 



where ttit and px are transverse mass and transverse 
momentum of the particle, respectively, and yb ~ 
\d.{^sn]si /ttin) is the beam rapidity [29| . 

According to the generalized blast-wave parametriza- 
tion, the radial flow rapidity which controls the magni- 
tude of the transverse expansion velocity is [gI. fiol Isol. . 



P = r{po + P2 cos(2(/)6)) 



(7) 



where r = y^{r cos{4>s)/ Rx)'^ + (r sm{(j)s)/ Ry)'^ ■ Po is 
the isotropic radial flow rapidity, and p2 is the amplitude 
of the anisotropic radial flow rapidity, respectively. The 
greater the magnitude of p2, the larger the momentum- 
space anisotropy. Here, (f>s is the azimuthal angle in co- 
ordinate space and ipb is the azimuthal angle of the boost 



r-ydyj- 
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In Fig. 2, the px spectra of six particles, n"^, K"^, p, 
and p, of the same sample, are presented by red solid cir- 
cles. The differential elliptic flow V2{pt) of pions, kaons, 
and protons are presented in Fig. 3 by black triangles, red 
solid circles and blue triangles, respectively. The error 
bars only include statistical errors. They are very small 
in comparison with the experimental data [l9| . Typically, 
the systematic errors are considered to be 5% when flt- 
ting the simulated data [131 • Due to resonance decays in 
the low momentum region of pions [20| . the data points 
in the low px regions of the spectra are excluded in this 
fltting. 

Using Eqs. (6), and (8), the fltting curves in each plots 
of Fig. 2 and 3 are drawn. They describe well the cor- 
responding data points of the px spectra and elliptic 
flow. The fltting parameters are T = 96.1 ± 1.0{MeV), 
po = 0.73 ± 0.01, and p2 = 0.035 ± 0.003. This temper- 
ature is the same magnitude as that given by hydrody- 
namics [^ [l3l , and experimental data [l^ • The param- 
eter of anisotropic radial flow rapidity, p2, is very close 
to the azimuthal amplitude of the suggested distribution, 
yT2 = 0.0334 ±0.0002. 
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FIG. 4: The centrality dependence of (a) the isotropic radial 
flow rapidity and (b) the temperature obtained by blast-wave 
parametrization from the sample of Au+Au collisions at 200 
GeV generated by the AMPT model with string melting. 



lated, as expected [l^, l35| . This is because the produced 
hadrons in central collisions have more time to cool down, 
rescatter, and form a stronger radial expansion flow. In 
peripheral collisions, there is not enough time to convert 
thermal energy to the collective flow motion. 



FIG. 2: (Color online) The transverse momentum spectra for 
TT^, , p and p within < 0.1 for the sample of Au+Au 
collisions at ^sjviv ~ 200 GeV generated by the AMPT model 
with string melting. 
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FIG. 3: (Color online) The differential elliptic flow V2{pt) 
for different particle species within < 0.1 for the sample 
of Au+Au collisions at ^/sNN ~ 200 GeV generated by the 
AMPT model with string melting. 

The centrality dependence of p2 is shown in Fig. 1(d) 
by black triangles. We can see that at each centrality, 
P2 is very close to yT2- The azimuthal amplitude of the 
suggested distribution coincides with the parameter of 
anisotropic radial flow rapidity extracted from a general- 
ized blast-wave parametrization. 

To complete the discussion of the fitting, the central- 
ity dependence of po and T are shown in Figs. 4(a) and 
(b), respectively. We can see that po increases with the 
increasing of number of participants, and reaches a maxi- 
mum in central collisions. The temperature changes with 
centrality in the opposite way as that of radial flow, po- 
The radial flow and the temperature are negatively corre- 



IV. SUMMARY AND CONCLUSIONS 

In this paper the measurement of the azimuthal distri- 
bution of mean transverse rapidity of final state hadrons 
is proposed as a more direct probe of the transverse mo- 
tion of the source than the known azimuthal multiplicity 
distribution. 

Using the sample generated by the AMPT model with 
string melting, we show that the isotropic part of the 
distribution is a combination of the radial and thermal 
motions. The azimuthally dependent part measures the 
anisotropy of transverse motion arising from non-central 
collisions. 

Using a generalized blast-wave parametrization, we 
further extract the temperature and radial flow param- 
eters from the same sample. It is found that the pa- 
rameter of anisotropic radial rapidity coincides with the 
azimuthal amplitude of the suggested distribution. 

The (yT(0)) provides a direct measurement of the 
anisotropic radial rapidity. This is important for hydro- 
dynamic calculations and for a direct measurement of 
shear viscosity in relativistic heavy ion collisions (36l | . 
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